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I. INTRODUCTION 



Multiple particle production in relativistic heavy-ion collisions is not fully understood yet. One of the most popular 
microscopic approaches to this phenomenon is formulated within the chromoelectric flux tube (or rather string) model 
. The flux tubes are assumed to be initially stretched between the constituents (quarks and diquarks) of the nucleons 
of colliding nuclei. As the constituents are flying away, the energy of a color tube increases making it unstable against 
the production of qq or dd pairs from the vacuum. New hadrons, created in the course of the color tube fragmentation, 
can also interact within hot and dense nuclear matter, and the in-medium cascade develops. 

The formation and subsequent break-up of color tubes/strings is a common feature adapted by all microscopic 
string models (see, e.g., and references therein), which are intended to describe dynamics of heavy- ion 

collisions at relativistic energies. For a uniform chromoelectric field E the probability to create a pair of quarks with 
mass TO, effective charge e, and transverse momentum pt per unit time and per unit volume reads 



P{pt)d^Pt = In |l - cxp 

The integrated probability is given by series expansion 
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which reproduces the classical Schwinger result Q derived in spinor quantum electrodynamics (QED) for e+e" 
production rate in the constant electric field. Usually, only the leading term in Eq. ^ is taken into account. 
According to this formula, the ratio of production rates of strange to nonstrange quark pairs, widely known as 
strangeness suppression factor 7^, is 
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where k = \eE\ is the so-called string tension. It appears that microscopic models underestimate the yield of 
strangeness in ultrarelativistic heavy-ion collisions j^. This issue is very important, because the abundant yield 
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of strange particles was predicted [3 as one of the signals of quark-gluon plasma (QGP) creation. According to 
Eq. the strangeness production can be enhanced either by taking into account strin g-string interaction, which 
leads to fusion of strings and formation of color ropes with larger effective string tension [T(l TTH fo ! IT^I l . or by dropping 
the quark masses 0|, e.g., due to chiral symmetry restoration. Also, the effects related to the finiteness of the strings 
can modify the production rates p"5t .16j . 

In the present Letter we explore another possibility: The real fields emerging in heavy-ion collisions act at the time 
interval comparable with the Compton scale. Since the dynamics of particle creation in time- dependent homogeneous 
fields differs from that of stationary fields [l7j . it is essential to properly modify the system description at early 
stage of nuclear collisions. The kinetic equation (KE) is a convenient tool to study the nonequilibrium evolution 
processes. The source term descri bing the vacuum pair creation process can be incorporated into the KE either 
in a phcnomcnological manner [TM Il9l | on the basis of Schwinger-like formula, or derived in a more sophisticated 
manner from the microscopic equations of motion pol | . Although it is believed that the phenomenological source term 
correctly reproduces qualitative features of the quantum mean field theory such approximation has yet to be 
verified for, e.g., time dependent fields or multicomponent systems. Compared to the phenomenological treatment, 
the approach within the framework of a transport equation "2^ contains several new dynamical aspects, such as 
longitudinal momentum dependence of the distribution functions and non-Markovian character of the time evolution. 
It takes into account effects of the field switching and statistics, as well. Therefore, the abundances of newly produced 
particles may considerably deviate from the values obtained for the constant and infinite field. The appearance of the 
non-Schwinger regime in a creation of electron-positron pairs in periodic laser field has been discussed in ^22, 23, 2^]. 
Another noteworthy feature of the modified source term is the suppression of zero momentum bosons |25|. This 
circumstance causes the "fermion dominance" effect at the short time scales j2^ and can lead to the abundant 
production of heavy fermions at expense of light bosons. 

In case of the QGP creation the characteristic time of the field variation is estimated to be of order of few fm/c 
|27j| and, therefore, the assumption of the space-time unvarying field is too crude. Hence, it is necessary to elaborate 
on the dynamics of particle vacuum production at short time scales, which are compatible with the particle inverse 
masses. Particularly, the dependence of the production rates on the shape and duration of the field pulse should be 
investigated. 

The paper is organized as follows. Kinetic equation with the source term, describing the vacuum production of 
fermions, is derived in the coUisionless limit in Sec. ^ Section imi presents study of the infiuence of the field pulse 
characteristics on the time-dependent vacuum production rates for the energies of RHIC and LHC. Comparison with 
the production rates, obtained for the infinite and constant chromoelectric field by the classical Schwinger-like formula, 
is also performed. Finally, conclusions are drawn in Sec. II VI 



II. KINETIC EQUATION 

As was mentioned above, the kinetic equation is a standard tool to describe the non-equilibrium evolution of a many 
body system. The general form of the KE for the distribution function /(p, t) in a strong spatially homogeneous time- 
dependent field is [13 

'-i^-,emf-^-Sip,t) + Cip,tl (4) 

where S{p,t) and C{p,t) are the source term and the collision integral, respectively. The source term describes the 
vacuum production of particle-antiparticle pairs in the external field, while the collision term governs their rescattering 
dynamics, which drives the system toward thermal equilibrium. For the sake of clarity, it is relevant to consider the 
coUisionless approximation, where the collision term in the r.h.s. of Eq. Q is omitted. 

The kinetic description of quarks in framework of the Wigner function formalism "2^, "2^ leads to a very complicated 
system of partial differential equations hard to solve both numerically and analytically. Therefore, following the 
approach of [TtL |2^ , we employ the canonical Bogoliubov transformation method within the QED approximation. To 
derive the source term for the system of fermions in an external electric field, we start from the Dirac equation 

{ij^'d^ - e-f^'Af, - m)i,{x) = . (5) 

Using the simple field configuration with vector potential in the Hamilton gauge = (0, 0, 0, A{t)) and homogeneous 
electric field E(i) — (0, 0, E{t)), E{t) = —A{t), one looks for the solutions of the Eq. in the form 

= K^o - 7k + 67^^(0 + m] x^^^k, t) Rr e^'^% (6) 
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where the superscript (±) denotes eigenstates with the positive and negative frequencies. Herein the spinors Rr 
(r = 1,2) are eigenvectors of the matrix 7*^7^ satisfying the condition R^Rs = ^S^s ■ The functions x^^H^^iO obey 
the oscihator-type equation 

xW(k,<)+K(k,t)-*ei(t)] x^^\k,t)^0, (7) 

where we define the total energy (k, i) = + (t) , the transverse energy = + , and the longitudinal 
momentum Puit) — fcy — eA{t). The solutions x*-^Hk, of Eq. Q for positive and negative frequencies are fixed by 
their asymptotic behavior at to = t ^ —00, where A{to) = 0. In this limit, the field operators ip{x) and ip{x) can be 
decomposed by the complete and orthonormalized set of spinor functions © as follows: 



4,{X) = J2 [^trH^) bM + ^iVi^) d\rito)] ■ (8) 



The operators &kr(io)i ^kr(^o) c^krC^o): c^kr(^o) describe the annihilation and creation of particles and antiparticles 
in the in-state |0i„) and obey the standard anticommutation rules. The time evolution leads to the mixing of states 
with positive and negative energies and, therefore, non-diagonal terms in the Hamiltonian corresponding to Eq. (jsj 
emerge. The diagonalization of the Hamiltonian, which is equivalent to the transition to quasiparticle representation, 
is performed by the time-dependent Bogoliubov transformation 

&kr(i) = auit) &kr(io) + /3k(t) d\rito) , 

with the imposed condition |ak(OI^ + |/3k(i)P ~ 1- The new operators 6kr(0 and dkr(i) describe the processes of 
quasiparticle creation and annihilation. By virtue of the Lagrange multipliers, one can find from the equations of 
motion ((7J) that the coefficients in the Bogoliubov transformation @ are connected via the relations ^3 



_ eE{t)e^ 2i9(k.to,t) 



where the dynamical phase is defined as 



(10) 



e(k,to,t) = / dt'uj{k,t') . (11) 

"'to 



To absorb the dynamical phase it is convenient to introduce new operators 

i?kr(t) = bi^rit) e-^»(k.*o.*) ^ 

D^t) = dk.(i) e-Wo,t)^ 



(12) 



which obey the anti-commutation relations: 

{Bkrit),B+^,it)} - {Dk,{t),D+^,{t)} = drr' <5kk' ■ (13) 

These operators satisfy the Heisenberg-type equations of motion 

dBk,.(<) _ eE{t)e^ 



dD.^rit) _ eE{t)e^ + 



(14) 



dt 2cj2(k,t) krV 

where H{t) is the Hamiltonian of the system of quasiparticles 

H{t) = ^(^' BUt) - D-^r{t) DtUt)] . (15) 

r,k 

The first term in the r.h.s. of Eqs. H14|l arises because of the unitary non-equivalence of the transition from the 
representation (jSJ to the quasiparticle one. 
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Next consider the evolution of the distribution function of quasiparticles with the momentum k and spin r defined 

as 

/.(k,t) = (0»„|5+ (t) 6krW|0„) = {Q^n\B+^.{t) Bkr{t)\0^n) • (16) 

According to the charge conservation the distribution functions for particles and anti-particles are related as /r(k, t) — 
/r(— k, t). Taking derivative in Eq. Hl()|) with respect to time t we have 

^ -Sit ^ ™ 

Here the function $r(k,i) = {Oin\D-'kr{t) i?kr(i)|Oi„) describes the vacuum production of pairs in external electric 
field E{t). Applying the equations of motion itT^ . one finds 

^ - lio) p^'<'''" - '1 - """'^'^ ^ <'^' 

The solution of Eq. ()18fl with the initial condition <i>r(k, to) = may be written in the following integral form 



t 

a>.(k, i) = Y / dt'^^^ [2/.(k, t') - 1] exp [2i0(k, t', t)] 

to 



(19) 



Inserting this result into the r.h.s. of Eq. p7|) we obtain the anticipated kinetic equation 

dfr{k,t) eE{t)e^ }^.eE{t')e, 



dt 2^2 (k^i) 



Since the distribution function does not depend on spin, the subscript r can be dropped: fr = /■ Substitution p = 
k — eA(t), where the 3-momentum is decomposed onto the transverse and longitudinal components p = p(pj_,P||(t)), 
yields to the reduction of the KE to Eq. with the source term 

5(p, t) = ^£;(t)w(p) j dh S(ti)w(p(i, tl)) [1 - 2/(p(t, tl), tl)] cos j dt2 Lo{t, h) j , (21) 

to \ tl / 

where w(p) = e^/a;(p) and 

^(t,ti) = ^el+p\{tM). 

p(t,ti) = p-e f E(t2)dt: 
Jti 



^2 • 



(22) 



The source term (I21f) demonstrates several interesting features, such as the dependence on particle longitudinal and 
transverse momentum, the account for spin and statistics, and the non-Markovian character of the time evolution. 
The me mory effects are caused by the time integration over the statistical factor (1 — 2/) and the non-local cosine 
function |25j| . while the structure of the coefficient w(p) defines the shape of the momentum distribution of created 
particles. 

In the coUisionless approximation the kinetic equation Q with the source term given by Eq. (|21|l can be transformed 
into the system of three ordinary differential equations 30J 

2/ — eEwvi , 

vi = e£:w(l - 2/) - 2ujV2 , (23) 
V2 = 2ljvi , 

where the dot denotes the full time derivative, and the auxiliary functions wi, V2 defined as 

vi = -2i?e$r(k,t), V2 = -2/TO<I>r(k,t). (24) 
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If the field strength is of order of the critical value, then it is necessary to take into account the back reaction of 
produced particles on primary field 0, |^ . The newly created particles polarize the vacuum and are accelerated by 
the external field. Their motion generates an internal field that in its turn modifies the initial background field. For the 
description of this phenomenon, the background field E{t) in l|21|l should be replaced by the sum E{t) = Eex{t)+Ei„ {t), 
where the generated internal field Ein can be found from the Maxwell equation |2Q|| 

^'" = -7?fe / — [2P„/+P.«i] , (25) 

with gf — 2Nc being the degeneracy factor and Nc the number of color degrees of freedom {Nc — 3). The total 
current density in the r.h.s. of Eq. H25() is the sum of conductivity and vacuum polarization currents, respectively. The 
integrand in (I25|l contains the logarithmic divergence which should be removed somehow by means of a regularization 
procedure 30]. We use here the simple ultraviolet cut-off of momentum integration on the border of a grid. The 
Eqs. (|23|l and H25|) represent the closed system of equations for the numerical analysis of the back reaction problem. 



III. VACUUM CREATION OF QUARKS WITH DIFFERENT MASSES 

The derived formalism can be applied to study the vacuum creation of quark- antiquark pairs in heavy ion collisions. 
Of special interest is the analysis of effects caused by the fast change of (chromo)electrical field at the time scale 
compatible with the inverse quark masses. To investigate a role of the field switching on/off effects, we approximate 
the time dependence of the flux-tube field by a short pulse oriented along the collision axis of primordial nuclei 

E{t)^E^c^p[-i2t/T)^] , (26) 

where Em is the field magnitude and r is the effective pulse width. The integer exponent n governs the steepness of 
the pulse: for n 1 (we restrict ourselves to Umax = 6), the pulse (|25|l becomes close to the rectangular one. The 
case with n — 2 reproduces well the soliton-like pulse 

E{t) = EmC0sh-\2t/T) , , . 

A{t) = -EmiT/2)tanh{2t/T) , ^ ' 

for which the Dirac equation allows exact analytical solution 0, • 

It is convenient to define the time dependent strangeness suppression factor fs as the ratio of densities of strange 
(s) to nonstrange (q) quarks 

, , nim^, t) 

Mt) = , 28 



where in the case of axial symmetry 



(2^ 



n{m-f,t)^gf / -—^f{mf,p^,p^^). (29) 



The system (|23|) is integrated by the Runge-Kutta method with the zero initial conditions /(p,to) — vi{p,to) — 
V2(jP-, to) = 0. The momentum dependence of the distribution function f{p\\,p±, t) is determined by coarse-graining of 
the momentum space to a 2-dimensional grid; in each of its node the system of equations (|23f) is solved. The parameters 
of the grid depend on the field strength, the typical values are Ap « 0.1m (step of the grid) and Pmax ~ 15 — 20to 
(boundary of the grid), so the total number of the solved equations is about 10®. 

Effect of the field pulse duration. To study the influence of the field pulse duration on the ratio of strange 
and nonstrange quark pairs we compare results, obtained for the time-dependent ratio H28|l with the soliton-like pulse 
(|27|1 . with those yielded for the constant chromoelectric field by Eq. ||2Jl. Input parameters, such as quark masses and 
string tension, are chosen according to |3|| , namely 

= 350 MeV , 
mu,d = 230 MeV , 

K = \eE\ = 0.9GeV/fm . (30) 

Since the critical field is defined as eEcr = rn^, the chosen value of the field is under-critical one for heavy strange 
quarks and over-critical one for the light nonstrange quarks. To solve the Maxwell equation H25|l the value of the 
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effective charge e should be determined. This can be done on basis of the estimates for the initial energy density e, 
which varies from 50 GeV/fm"^ at RHIC to 520 GeV/fm'^ at LHC [s^l- Assuming for simplicity that all initial energy 
is deposited in the field sector, one obtains E^x ~ 3.16 GeV/fm for RHIC and E^x ~ 10 GeV/fm for LHC. The 
corresponding values of the effective charge are e w 0.22 and e w 0.07 for RHIC and LHC, respectively. 

The results of calculations of the asymptotic value fs{t ^ r) of the suppression factor fs are depicted in Fig. ^ for 
the field pulse (|27|l . Trivially, the Schwinger-like estimate for the constant and infinite field © gives the constant 
value fs = 0.3. The time-dependent case is more complex: for pulses shorter than the particle inverse masses the 
creation of ss quark pairs is significantly enhanced. This is a direct consequence of the uncertainty relation for the 
energy and time. It is worth mentioning that already at r > 3 fm/c the exact result is very close to that given by 
the Schwinger formula ©■ However, for the characteristic time of QGP formation Tq ~ 1 fm/c js^, the production 
probability of strange quarks in the soliton-like field yield (|27|l is at least 1.5 times larger than that in the stationary 
case. 

Influence of the field pulse shape. Figure |2 shows the time evolution of the factor fs{t) during action of field 
pulses with n — 6 and f?7| with the same width r = 1 fm/c. We see, that production of strange quarks increases 
with rise of the power n in exponent in Ea. (|26|l . Note also, that the intermediate effective value of the time dependent 
ratio fs is considerably larger than its final one: it varies from about 0.6 to 0.39 for the soliton-like pulse and from 
0.75 to 0.52 for the rectangular one. The evolution of densities of the created particles is presented in Fig. |21 The 
density of both light and heavy quark pairs in the pulsing field initially increases. As the field saturates and starts to 
decrease, the process of particle absorption by the field dominates the particle production one. The longer the field 
pulses, the stronger absorption. Therefore, the way of the field oscillation can significantly change the value of the 
suppression factor. The similar result concerning the role of pulse shape of periodic laser field on the electron-positron 
vacuum production rate was obtained in |23j within the framework of the approximate imaginary time method. 

In Ref. 13j two ways of increasing the strangeness production within the framework of the Schwinger mechanism 
have been discussed, namely (i) either the field strength (string tension) is increasing, or, equivalently, (ii) quark 
masses are dropping due to chiral symmetry restoration. Assuming a significant reduction of quark masses from H^Ojl 
to the current quark masses rus = 0.23 GeV, iriq = 0.01 GeV 13], we get from Eq. J^Jl the enhancement of strangeness 
production in a constant infinite field from 7^ « 0.3 to 7^ w 0.4, which is equivalent to increase of the string tension 
K from 0.9 GcV/fm to 1.22 GeV/fm. 

Let us perform similar fitting procedure for the rectangular pulse given by Eq. H26|) . The corresponding dynamical 
picture is displayed in Fig. ^ It appears that the scenarios considered above are not fully equivalent in dynamical 
sense, because the production of light particles at very short times is apparently suppressed. Therefore, the reduction 
of quark masses due to the chiral symmetry restoration is even more effective for the enhancement of strangeness 
production than the mechanism of string overlap, which leads to the formation of color rope and increase of the string 
tension. 

The solution of Eas. (|23|) and H25|l shows that the back reactions play minor role in the pair production processes 
with the values of parameters given by the set (|^ . because the current of secondaries is quite weak. The increase 
of the string tension k from 0.9 GeV/fm to 1.22 GeV/fm or the reduction of quark masses to nis — 0.23 GeV and 
rUq = 0.01 GeV does not significantly modify this scenario. 

IV. CONCLUSIONS 

Kinetic equation with the source term, describing the vacuum production of fermions in a strong, time-dependent 
field, is used for the two-component system investigation in the coUisionless approximation. It is shown that dynamics 
of vacuum particle creation at time scale compatible with the inverse mass of particle depends essentially on the specific 
characteristics of field configuration (a form and a duration of the field pulse). In particular, the Schwinger-like regime 
of particle creation might not be realized at the typical time scale of QGP formation, tq 1 fm/c. As a consequence, 
production of heavy strange quarks becomes more abundant. The role of string fusion and reduction of quark masses 
due to the chiral symmetry restoration which can alter the strangeness production is studied. Within the proposed 
dynamical scenario the mass reduction mechanism appears to be more effective for the enhancement of strangeness 
yield than the formation of color rope. Since the current of the produced secondaries is weak, the contribution of the 
back reactions to the production process is small. Our study shows that the time evolution picture of the gluon field 
should be incorporated consistently into the color flux tube model for quantitative description of parton production 
in ultrarelativistic heavy ion collisions. 
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FIG. 1: Pulse width dependence of the asymptotic value fs{t S> t) of the suppression factor fs, the dotted hue shows the 
corresponding estimation from the Schwinger formula Q. 
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FIG. 2: Time dependence of the suppression factor fs{t) for two types of field pulse, given by Eq. I|26|l . witli the same width; 
the pulse with n — 6 (solid line) and the pulse with n — 2 (dash-dotted line). Dotted and dashed lines show the corresponding 
field profiles. 
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FIG. 3: Time dependence of particle densities normalized on its final value Ur = n{t ^ r) for pulse 1271 with the width to = 1 
fm/c. 




FIG. 4: Comparison of two mechanisms of fs enhancement: the quark mass reduction (sohd hne) and the string tension 
increase (dotted hne). 



